Fermi surface shrinking and interband coupling in iron-based pnictides 
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Recent measurements of Fermi surface with de Haas-van Alphen oscillations in LaFePO showed 
a shrinking of the Fermi pockets with respect to first-principle calculations, suggesting an energy 
shift of the hole and electrons bands with respect to local density approximation. We show that 
these shifts are a natural consequence of the strong particle-hole asymmetry of electronic bands 
in pnictides, and that they provide an indirect experimental evidence of a dominant interband 
scattering in these systems. 
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A new challenge in the field of condensed matter is 
represented by the recent discovery of high-T c supercon- 
ductivity in the iron-based pnictide family First- 
principle calculations have soon identified four main 
bands at the Fermi level with two-dimensional charac- 
ter: two hole-like pockets around the T point and two 
electron-like pockets around the M point 0, 0, 0, H, S 0] • 
A fifth three-dimensional band crossing the Fermi level 
at the r point is sometimes predicted by Local Density 
(LDA) calculations @, 0], even if its position with re- 
spect to the chemical potential is strongly affected by 
the interlayer distance. Since determining the Fermi sur- 
face topology is the first step toward the understanding 
of these new materials, a lot of experimental work has 
been devoted to its investigation. A momentum-resolved 
mapping of the dispersion of the occupied quasi-particlc 
states in the normal state and in the superconducting 
one has been provided by angle-resolved photocmission 
spectrosc^y^ARPES), both in 1111 j| and 122 sys- 
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ing the Fermi surface topology, which is not momentum 
resolved but has the advantage of being a bulk probe, 
is based on dc Haas-van Alphen (dHvA) magnetization 
measurements, which allow one to estimate the size of 
the Fermi areas and the effective mass m* for each Fermi 
sheet 14, ]1| 16, 17, H)]. Although an overall qualitative 
agreement is found between LDA calculations and these 
experiments, interesting discrepancies still remain. In 
general, we can distinguish between high energy and low 
energy discrepancies. ARPES data in Bao.6Ko.4Fe2AS2 
El and LaFePO [8] for instance find overall band- 
widths a factor 2 smaller than the LDA calculations, pos- 
sibly related to the static electronic correlation [jjj . On 
the other hands, the ARPES observation of Fermi veloci- 
ties a factor 4 larger than the LDA calculated low-energy 
ones [10| points out a relevant role of dynamical quan- 
tum renormalization effects associated with a low-energy 
boson-mediated interaction. 

Besides the bandwidth renormalization, a second strik- 
ing result emerging from dHvA is a substantial shift of 
the bands with respect to the Fermi level when compared 



to LDA calculations. Indeed, the accurate determination 
of the Fermi surface areas provided by dHvA gives values 
smaller then those expected by LDA. Such discrepancy 
can be accounted for assuming a shift of the LDA bands. 
Notably, such shifts have different signs in hole and elec- 
tron bands, being downward for the hole-band and up- 
ward for the electron ones 3, HI • This effect seems to 
persist to higher energies, as one can infer from ARPES, 
which gives a systematical reduction of the energy differ- 
ence between the bottom of the electron bands and the 
top of the hole bands with respect to the values predicted 
by LDA 0, EH, EH, EIJ ■ Even though a certain degree 
of inaccuracy could be present in density-functional the- 
ory calculations, the persistent observations of such shifts 
suggests a robust feature in these materials whose origin 
is at the moment unknown. 

In this Letter we present a comprehensive explanation 
for the origin of the band shifts observed in dHvA ex- 
periments. We show that they are a natural consequence 
of the multiband character of these systems and of the 
strong particle-hole asymmetry of the bands, which in- 
duce finite-band self-energy effects that are usually irrel- 
evant for the half- filled single-band case. We also show 
that a simple analysis of the dHvA data provides evi- 
dence that the dominant interaction in these systems is 
the interband one between hole and electronic states, as 
it has been already argued on the basis of the nestin g be - 
tween hole and electron Fermi sheets @, 0, HJ [53^23] • 
More precisely, our calculations give an estimate of the 
interband coupling V of the order V « 0.46 eV, which 
gives rise also to a mass enhancement Z a w 1.4 for each 
band. 

We start our analysis by discussing first a generic 
multiband system, where the electrons are coupled 
to a bosonic mode with local propagator D(lji) — 
f (K}2flB(ri)/(ri 2 +ujf ), and where B(Q) is the density 
of states of the bosonic excitations. The self-energy in 
the Matsubara space can be written as: 



Y, a (iu n ) = -TyV at pD{u) n - u m )Gp{iu m ), (1) 
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where a, (3 are band indexes, V a ,p (V a ,p = Vg, Q ) is the 
multiband interaction, which for a retarded interaction is 
always positive (V a .p > 0), and G a (z) is the local Green's 
function for the a band, namely: 



G a (z) 



deN a (e) 
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z — e — T, a (z) + [i 
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Here N a (e) is the electronic density of states of the a 



band, with upper and lower band edges E n 



and 



E m i n , a , respectively, and we drop the spin index since 
it does not play any role in the following. 



In the conventional Eliashberg analysis [24J, one usu- 
ally assumes that the distance of the chemical potential 
from the band edges is much larger than the typical bo- 
son energy scale, so that one can approximate the DOS 
with its value at the Fermi level, N a (e) N a (fj,), and 
one can extend the integration limits over e in Eq. |2]) to 
±oo. In this way one is enforcing the particle-hole sym- 
metry and the Matsubara self-energy T, a (iuj n ) is purely 
imaginary. Enforcing implicitly the particle-hole symme- 
try can be however quite dangerous in systems as the 
iron-pnictides where each band is strongly away from 
the half-filling. As a byproduct of taking into account 
the particle-hole asymmetry, the self-energy acquires a 
finite real part Xa(iu n ) = R.e£ Q (icj„) ^ 0, whose low 
energy limit \ a — Xa{i^n=o) gives rise to a band shift 
that can be in general different in each band. While this 
effect is usually disregarded in single-band systems, be- 
cause it can be absorbed in a redefinition of the chemical 
potential, in the multiband case it can lead to observable 
relative shifts of the various bands with respect to the 
Fermi level. 

Before undertaking a fully numerical solution of Eqs. 
((T|)-([2|), it is instructive to consider the result of the lowest 
order perturbation theory, where the Green's function 
Gf3(iu> m ) in Eq. |T]) is taken to be the non- interacting 
one with S ( g(io; m ) = 0. We shall assume also purely two- 
dimensional parabolic bands, so that N a — l/(-E m ax,a — 
£"min,a) and we consider for the moment, for simplicity, 
a Einstein mode B(Q) = (ojq/2)S(Q — ujq). We get thus 
the analytical expression (for T « 0): 
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In most cases, the exchanged boson energy is the lowest 
energy scale in the system, so that 



Xa ~ -^y^Va^Npln 
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Eq. ((!]) gives correctly \a = in the case of particle-hole 
symmetry \E m&Xi 

Ml — l-^min On the other hand, 

in an electron-like band \E max ,f3 — fi\ > \E m i n ^ — while 
the opposite inequality applies for a hole-like band. Let 
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FIG. 1: (Color online). Sketch of our band structure model 
(a) and of the corresponding Fermi surfaces (b). We also show 
the only non-zero interband couplings. 



us consider now a single band system. In this case the 
band shift \a,o is compelled to be positive (upward) for 
a hole-system and negative (downward) for an electron- 
like band. Such observation seems at odds with what 
reported in LaFePO. 

Things are however more subtle in multiband systems, 
where the shift of each band depends also on the con- 
tribution of all the other bands weighted with the corre- 
sponding interband coupling. This effect is particularly 
important in pnictides where the dominant coupling is 
thought to be the spin-mediated interband interaction. 
In such a situation, the particle-hole asymmetry of the 
electron bands is responsible for the downward shift of 
the hole-like bands which, vice versa, give rise to the up- 
ward shift of the electron ones. We would like to stress 
that these argumentations do not depend on the band 
or interaction details. In full generality the simple ob- 
servation of a upward shift of the electron bands and the 
downward shift of the hole-like ones is a direct experi- 
mental evidence in these compounds of the dominance of 
the interband coupling on the intraband one, which would 
produce the opposite scenario. 

Within this framework, a more detailed analysis of 
the dHvA 



measurements in LaFePO 14| can provide 
an useful insight into the strength of the interband cou- 
pling in pnictides. Following Ref. [25} we consider four 
bands (1,2 hole-like and 3,4 electron-like) and we as- 
sume, because of the nesting properties, a purely in- 
terband scattering connecting hole and electron bands, 
namely Vj.,3, V1.4, V 2 ,3, V 3 ^, (see Fig. [T]). Note that, un- 
like to other pnictides, all the Fermi areas in LaFePO 
are quite comparable, so that we can expect that in- 
terband scattering V a ^ does not depend much on the 
band index, and we can assume, as a first approxima- 
tion, Vi.,3 = Vi,4 = V<2,3 = V34 = V. We model each band 
with a purely two-dimensional parabolic dispersion, and 
we assume for simplicity bands 3, 4 to be completely de- 
generate. The band parameters are extracted from the 
LDA calculations 0. More explicitly, focusing on the 
kz = plane, we determine the Fermi vectors k-p, a from 
the T-X cut for the hole-bands, and from the M-X cut for 
the electron ones. Setting for convention the Fermi level 
fi = 0, we also take from LDA calculations the nearest 
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LDA calculations LaFePO 




1 


0.58 0.205 -5.031 0.191 
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1.0 
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1.14 0.205 -2.462 0.375 
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0.79 3.551 -0.295 0.260 
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Renormalized LDA calculations 


LaFePO 




1 


1.16 0.102 -2.516 0.382 
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2.28 0.102 -1.231 0.750 


-26.5 


3.2 


3,4 


1.58 1.776 -0.147 0.520 


31 
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TABLE I: Microscopic band parameters extracted from un- 
renormalized LDA calculations 2] (top), and renormalized 
ones (bottom). Also shown are the calculated band shifts 
Xo,a and the renormalized mass m* for V = 1.55 eV and 
V = 0.46 eV, respectively. m e is the free electron mass. 



band edge for each band, namely £7 m ax,i> £max,2 for the 
hole-like bands, -Emm, 3, -Emm, 4 for the electron ones, and 
we use it to estimate the non-interacting mass for each 
band m a = h 2 k^ ct /2£ , max ( min ) a and the corresponding 
density of states N a = m a a /2tt?i 2 ( a is here the in- 
plane lattice constant). Finally, the effective band edge 
far from the Fermi level [26j follows from the relation 
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V (eV) 
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V (eV) 



N, 



l/(i? maXiQ — i?min,a)- All the estimated values are 
reported in Tabic [U 

We can now employ a full numerical solution of Eqs. 
([I])-© to get a qualitative estimate the strength of the 
interband coupling V by comparing the calculated mag- 
nitude of the band shifts and the effective mass to* 
with the experiments. This latter quantity can be eval- 
uated from the numerical solution of Eqs. (P)-© as 
to* = m a Z a , where Z a = 1 — lmY* a (iu> n ) /Lu n \ n= o- To ac- 
count for a spin-mediated interaction mechanism we use 
the Lorentzian spectrum typical of spin fluctuations 27 1 
B(fl) oc £1(jJo/it(oJq + f^), with the characteristic energy 
scale uj = 20 meVjH HE [3l|. The smallness of 
ujo compared with the Fermi energies, E-p 1=3 100 — 200 
meV (Table H|, guarantees the validity of Migdals' the- 
orem 1321 also for spin- fluctuations, as discussed e.g. in 

Refs. [H IH 13. In Fi S- I3 1 we show the band shifts 
Xa,o evaluated from the numerical solution of Eqs. |(TJ)- 
© as a function of V by using the set of microscopic 
parameters obtained by bare LDA calculations. From 
Fig. [2^, we get an estimate V ~ 1.55 eV to account 
for the hole band shift A12 ~ 53 meV reported in Ref. 
14j . This value V ~ 1.55 eV gives also a theoreti- 
cal positive shift X3 = X4 ~ 60 meV for the electron 
bands, which agrees quite well with the average exper- 
imental value A 3j4 = (A 3 + A 4 )/2 = 56.5 meVQ. It 
is clear that possible differences between the set of hole 
and electron bands can arise in a more refined treat- 
ment when the two electron bands are no more assumed 
to be degenerate, and when interband coupling is not 
taken constant. This value V ~ 1.55 eV would also 



FIG. 2: (Color online), (a) Bands shift Xo,a as function of 
the interband coupling V for the parameter set taken from 
LDA calculations (table U top). The horizontal dashed lines 
mark the average band shift as estimated in Ref. fl4l . fl5| ; 
(b) same as in panel (a) but considering LDA band structure 
renormalized by a factor 2 (table U bottom). 



give weak-coupling multiband dimensionless interactions 
A a ,/3 = Va^Np: A13 = A23 = A14 = A 2 4 = 0.40, 
A31 = A41 = 0.30, A32 = A42 = 0.58 which yield effective 
masses m *„lr n e . ~ 1, 1.4, 2 smaller than what reported in 
Refs. IHfillS. A numerical solution of the multiband 
Eliashberg theory would give also a critical temperature 
T c w 29 K, larger than the experimental one, T c ss 6 K. 

We identify the origin of this discrepancy in the lack, in 
the above analysis, of high-energy band renormalization 
pointed out by ARPES experiments [H, which show that 
the real band structure is roughly twice narrower than 
the LDA one. The nature of such high-energy renormal- 
ization is at the moment unknown, but it must be of 
different origin than the interband spin-fluctuation cou- 
pling which acts on the energy scale Uq. One can easily 
see indeed that for a retarded interaction the band renor- 
malization function is negligible Z{ui) ~ 1 when w>wo 
[24] ] . This means that the responsible mechanism for such 
high-energy renormalization must be operative on an en- 
ergy scale larger than the electronic bandwidth, suggest- 
ing once more a correlation-driven origin [19( . Note also 
that in this case such mechanism is not expected to affect 
the band shift since, as shown in Eq. ((4]), Xa,o w for 



UJQ > E, 



max(min) ,a 



We take into account the bandwidth narrowing by 
renormalizing the whole LDA band structure as ex- 
perimentally pointed out by ARPES in LaFePO and 
Bao.6Ko.4Fe2AS2- The corresponding microscopic pa- 
rameters in our model are listed in Table HI Note that, in 
such rescaled band structure, also the band shifts needed 
to recover the experimental Fermi areas measured by 
dHvA techniques are reduced by a corresponding factor 
2. The plot of the band shifts as function of the interband 
coupling, for such renormalized LDA band structure, is 
shown in Fig. Our estimate of the interband in- 
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FIG. 3: (Color online). Intensity map of the spectral function 
for the interacting hole bands (left panel) and the two degen- 
erate electron bands (right panel), obtained by using LDA 
parameters renormalized by a factor 2 and V = 0.46 eV. The 
dashed lines represent the non-interacting parabolic bands, 
and the horizontal solid line is the chemical potential. The 
ticks on the top mark the value of the Fermi vectors in the 
presence (solid yellow ticks) and in the absence (white dashed 
ticks) of interaction. Note the shrinking of the Fermi surfaces 
due to the coupling to the retarded bosonic mode. 



teraction in this case is thus V « 0.46 eV, which gives 
weak-coupling values A13 = A23 = A14 = A24 = 0.24, 
A31 = A41 = 0.18, A 32 = A 42 = 0.34, with a T c w 9 K 
and to*. / m e w 1.6 — 3.2, in better agreement with dHvA 



and specific heat measurements [15J, |35|, [3 



The corresponding intensity map of the spectral function 
for the interacting hole and electron bands, as obtained 
by the Marsiglio-Schossmann-Carbotte analytical contin- 
uation 37], is also shown in Fig. [3] along with the bare 
band dispersions renormalized by the factor 2. Here it 
is clearly visible the shrinking of the Fermi area of each 
band due to the interband self-energy effects (xa) which 
give rise to a downward shift of the hole bands and to 
an upward shift of the electron ones. Note that such 
self-energy corrections due to the particle-hole asymme- 
try survive until energies much larger than uiq, so that 
the energy shift is effective also at momenta far from 
the Fermi level. In particular this means that, when 
interband interactions are predominant, the top of the 
hole bands (at T) and the bottom of the electron bands 
(at M) approach each other with respect to the predic- 
tion of LDA, in agreement with the ARPES observation 
in 122 compounds [To|, [HL E3, Note also in Fig. [3] 

the effective mass renormalization at the Fermi velocity 
to* = m a Z a due to the retarded interband interaction. 
This effect, unlike the band shift Xa, disappears at the 
scale energy lvq, giving rise to a "kink" in the electronic 



dispersion, recently observed in pnictides 12j and widely 
discussed in the past in the context of cuprates [3^ . 

In conclusion, we demonstrated that the band shifts 
reported in pnictides when comparing the experimen- 
tally measured Fermi surfaces and band dispersions with 



LDA calculations are a direct consequence of the cou- 
pling to a bosonic mode, once that the strong particle- 
hole asymmetry and the multiband character of these 
systems are properly taken into account. Moreover, we 
showed that the sign of the measured shifts, with both 
the hole and electron bands approaching the Fermi level, 
suggests that interband interactions dominate over in- 
terband ones. Our results pose strong constraints on 
the theoretical modeling of interactions in pnictides, and 
show unambiguously that finite-band effects cannot be 
disregarded in these systems. 

This work was partially funded by MIUR project PRIN 
2007FW3MJX. 
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